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ABSTRACT 



> 



Aims. We present a detailed analysis of a very high resolution (R « 112, 000) spectrum of the quasar HE 0515-4414 obtained using the High 
Accuracy Radial velocity Planet Searcher (HARPS) mounted on the ESO 3.6 m telescope at the La Silla observatory. The main aim is to use 
HARPS spectrum of very high wavelength calibration accuracy (better than 1 mA), to constrain the variation of a = e 2 /he and investigate 
I any possible systematic inaccuracies in the wavelength calibration of the UV Echelle Spectrograph (UVES) mounted on the ESO Very Large 
■ Telescope (VLT). 

\Q Methods. A cross-correlation analysis between the Th-Ar lamp spectra obtained with HARPS and UVES is carried out to detect any possible 
. shift between the two spectra. Absolute wavelength calibration accuracies, and how that translate to the uncertainties in Aa/a are computed 
using Gaussian fits for both lamp spectra. The value of Aa/a at z a b s = 1.1508 is obtained using Many Multiplet method, and simultaneous 
Voigt profile fits of HARPS and UVES spectra. 

Results. We find the shift between the HARPS and UVES spectra has mean around zero with a dispersion of a = 1 mA. This is shown to 
O be well within the wavelength calibration accuracy of UVES (i.e a <a 4 mA). We show that the uncertainties in the wavelength calibration 
induce an error of about, Aa/a < 10~ 6 , in the determination of the variation of the fine-structure constant. Thus, the results of non-evolving 



Or 



ffi ' Aa/a reported in the literature based on UVES/VLT data should not be heavily influenced by problems related to wavelength calibration 
. . uncertainties. Our higher resolution spectrum of the z a t>s = 1 1508 Damped Lyman-fy system toward HE 0515-4414 reveals more components 
compared to the UVES spectrum. Using only Fe n lines of z abs = 1.1508 system, we obtain Aa/a = (0.05 ± 0.24) x 10~ 5 . This result is 
consistent with the earlier measurement for this system using the UVES spectrum alone. 
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1. Introduction 

Some of the modern theories of fundamental physics, such 
as SUSY, GUT and Super-string theory, allow possible space 
and time variations of the fundamental constants, thus moti- 
vating an experimental search for such a variation (Uzan 2003 
and 2004 for a detail review on the subject). Murphy et al. 
(2003), applying the Many Multiplet method (MM method) 
to 143 complex metal line systems, claimed a non-zero vari- 
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ation of the fine-structure constant, a = e 2 /flc: (Aa/a) 
( 0.57 + 0.11) x 10~ 5 for 0.2 < z < 3.5, where Aa/a = 
(a z - ao)/ao, with ao being the present value and a z its value 
at redshift z. This result, if true, would have very important im- 
plications to our understanding of fundamental physics and has 
therefore motivated new activities in the field. Search for the 
possible time-variation of a using alkali doublets has started 
long ago (Bahcall et al. 1967). The alkali-doublet method is a 
clean method for constraining the variation in a using spectral 
lines because it uses transitions from the same species (Wolfe 
et al. 1976; Levshakov 1994; Potekhin et al. 1994; Cowie & 
Songaila, 1995; Varshalovich et al. 1996; Varshalovich et al. 
2000; Murphy et al. 2001a; Martinez et al. 2003; Chand et al. 
2005). The tightest constraint obtained using this method till 
date is Aa/a = (0.15±0.44) xl0~ 5 at z ~ 2 (Chand et al. 2005). 
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Studies based on heavy element molecular absorption lines 
seen in the radio/mm wavelength range are more sensitive than 
that based on optical/UV absorption lines. They usually pro- 
vide constraints on the variation of a combination of the fine- 
structure constant, the proton g-factor (G p ) and the electron- 
to-proton mass ratio QS). Murphy et al. (2001b) have obtained 
Aa/a = (-0.10 + 0.22) x 10~ 5 at z = 0.2467 and Aa/a = 
(-0.08 ± 0.27) x 10~ 5 at z = 0.6847, assuming a constant pro- 
ton g-factor (G p ). It has been pointed out that OH lines are 
very useful in simultaneously constraining various fundamen- 
tal constants (Chengalur & Kanekar 2003; Kanekar et al. 2004; 
Darling 2003, 2004). These studies have provided Aa/a = 
(0.6 ± 1.0) x 10~ 5 for an absorption system at z a bs= 0.247 to- 
ward PKS 1413+135. Such studies have not been performed 
yet at higher redshift (i.e z > 1) due to the lack of molecular 
absorption systems. 

Constraints on the variations of a are also obtained from 
terrestrial measurements. The most stringent constrain has been 
obtained from the analysis of the Oklo phenomenon. Fujii et 
al. (2000) find that Aa/a = (-0.8 ± 1.0) x 10~ 8 over a pe- 
riod of about 2 billion years (or z — 0.45). Laboratory exper- 
iments also give very stringent constraints on the local vari- 
ation of a. Marion et al. (2003) have obtained, Aa/aAt = 
(-0.4 ± 16) x 10 _16 yr -1 , by comparing the hyperfine transi- 
tion in 87 Rb and 133 Cs over a period of 4 years assuming no 
variation in the magnetic moments. Fischer et al. (2004) have 
obtained, Aa/aAt = (-0.9+2.9)x 10~ 16 yr~', by comparing the 
absolute 15-25 transition of atomic hydrogen to the ground 
state of Cesium. A linear extrapolation gives a constraint of 
-1 .3 x 10~ 6 < Aa/a < 1 .9 x 10~ 6 at z = 1 for the most favored 
cosmology (£l m = 0.27, Q A = 0.73 and h = 0.71). 

Clearly all the experimental results summarized above are 
consistent with no variation of a. However, these results do 
not directly conflict with the positive detection by Murphy et 
al. (2003) either because of the insufficient sensitivity of the 
method (as in the case of alkali doublets) or because of the 
different redshift coverage (as in the case of radio and terres- 
trial measurements). However, recent attempts using the MM 
method (or its modified version) applied to very high quality 
UVES spectra have resulted in null detections. The analysis 
of Fe ii multiplets and Mg n doublets in a homogeneous sam- 
ple of 23 systems has yielded a stringent constraint, Aa/a = 
(-0.06 ± 0.06) x lO -5 (Chand et al. 2004; Srianand et al. 2004). 
Modified MM method analysis of z abs = 1.1508 toward HE 
0515-4414 that avoids possible complications due to isotopic 
abundances has resulted in Aa/a = (0.01 + 0.17) X 10~ 5 (Quast 
et al. 2004). Levshakov et al. (2005b) have re-analysis this sys- 
tem using the single ion differential alpha measurement method 
as described in Levshakov et al. (2005a), and obtained Aa/a = 
(-0.007±0.084)xl0- 5 . Clearly all studies based on VLT-UVES 
data are in contradiction with the conclusions of Murphy et al. 
(2003). 

A first possible concern about these studies is the real ac- 
curacy and robustness of the various calibration procedures. A 
second possible source of uncertainty comes from the multi- 
component Voigt-profile decomposition. It is very important to 
check how sensitive the derived constraints are to the profile de- 
composition. This can be done by performing the analysis on 



data of higher resolution than typical UVES (or HIRES) spec- 
tra. The best way to investigate all this is to compare data taken 
by UVES (or HIRES) with data on the same object taken with 
another completely independent, well controlled, and higher 
spectral resolution instrument. The advent of HARPS mounted 
on the ESO 3.6 m telescope makes this possible. Unfortunately 
this is only possible on the brightest quasar in the southern sky, 
HE 0515-4414. 

This forms the basic motivations of this work. We report 
the analysis of the z a bs= 115 DLA system toward QSO HE 
0515-4414 (De la Varga et al. 2000, Quast et al. 2004, 2005) 
using very high resolution (R~ 112, 000) spectra obtained with 
HARPS mounted on the ESO 3.6 m telescope. The organiza- 
tion of the paper is as follows. The HARPS observations of 
HE 0515-4414 are described in Section 2. Calibration accu- 
racy and comparison with the UVES observations are discussed 
in Section 3. In Section 4 we present the joint analysis of the 
HARPS and UVES spectra. Results are summarized and dis- 
cussed in Section 5. 

2. Observations 

The spectrum of HE 0515-4414 used in this work was ob- 
tained with the High Accuracy Radial velocity Planet Searcher 
(HARPS) mounted on the ESO 3.6 m telescope at the La Silla 
observatory. HARPS is a fiber-fed spectrograph and is there- 
fore less affected by any fluctuation in the seeing conditions 
(Mosser et al. 2004). It is installed in the Coude room of the 
3.6 m telescope building and is enclosed in a box in which vac- 
uum and constant temperature are maintained. The instrument 
has been specifically designed to guarantee stability and high- 
accuracy wavelength calibration. 

The observations were carried over four nights in classical 
fiber spectroscopy mode, with one fiber on the target and the 
other on the sky. The CCD was read in normal low readout 
mode without binning. The echelle order extraction from the 
raw data frame is done using the HARPS reduction pipeline. 
The error spectrum is computed by modeling the photon noise 
with a Poisson distribution and CCD readout noise with a 
Gaussian distribution. The calibrated spectrum is converted to 
vacuum wavelengths according to Edlen (1966) and the helio- 
centric velocity correction is done manually using the dedi- 
cated MIDAS (ESO-Munich Image Data Analysis Software) 
procedure. Special attention was given while merging the or- 
ders. While combining overlapping regions, higher weights 
were assigned to the wavelength ranges toward the center of 
the order compared to the one at the edges. The resulting 1-D 
spectrum covers the wavelength range from 3800 to 6900 A, 
with a gap between 5300 to 5330 A caused by the transition 
between the two CCDs used in HARPS. In total, we obtained 
14 individual exposures, each of duration between 1 and 1.5 
hour. Combination of individual exposures is performed using 
a sliding window and weighting the signal by the errors in each 
pixel. The final error spectrum was obtained by adding quadrat - 
ically in each pixel the extracted errors and the rms of the 14 
individual measurements. The final combined spectrum has a 
S/N ratio of about 30 to 40 per pixel of size ~0.015 A and a 
spectral resolution of Rx 112,000. 
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To make quantitative comparisons, as will be discussed in 
the next section, we have also used the UVES spectrum of this 
QSO. The details of the UVES observation and data reduction 
can be found in Quast et al. (2004). However we have used 
our procedures for air-to-vacuum wavelength conversion, he- 
liocentric velocity correction and for the addition of individual 
exposures as in the case of the HARPS spectrum. 

3. Accuracy of wavelength calibration 

In this Section we investigate (i) the cross-correlation between 
the Th-Ar lamp spectra obtained with HARPS and UVES, (ii) 
the absolute wavelength calibration accuracies of HARPS and 
UVES and (iii) how the uncertainties in the wavelength cali- 
bration translate into uncertainties in Aa/a measurements in 
the case of HARPS and UVES. 

3.1. Cross-correlation of UVES and HARPS Th-Ar 
spectra 

To estimate how well the UVES and HARPS wavelength scales 
agree, one can in principle use the narrow heavy element ab- 
sorption lines seen in the spectra of the QSO. However not 
only the number of such lines is small but also, due to differ- 
ences in the resolutions and S/N ratios, spurious shifts can be 
introduced in the analysis. In order to avoid this, we perform a 
cross-correlation analysis between the Th-Ar lamp spectra ob- 
tained with UVES and HARPS. We have 4 and 14 Th-Ar lamp 
exposures respectively for UVES and HARPS observations in 
the setting that covers the wavelength range where Fe n and 
Mg ii absorption lines from the z a b. s = 1.1508 absorption sys- 
tem are seen. We have combined all the extracted Th-Ar ex- 
posures after subtracting a smooth continuum corresponding to 
the background light. 

The cross-correlation analysis was performed on groups of 
five consecutive unblended emission lines that are clearly seen 
in both the UVES and HARPS spectra. For this, both spectra 
were re-sampled to an uniform wavelength scale using cubic 
spline and the pixel-by-pixel cross correlation was performed 
by shifting the UVES spectrum with respect to the HARPS 
spectrum. The results of the cross-correlation at places where 
absorption lines at z a b s = 1.1508 are redshifted are shown in 
Fig. n All the curves shown in this figure have their peak at 
zero pixel shift with a typical pixel size of 15 mA. In order to 
derive sub-pixel accuracy in the cross-correlation, we have fit- 
ted a Gaussian to the cross correlation curves as is shown by 
dotted lines (Fig.^ an d derive its centroid accurately. The cor- 
responding values are given in each panel. The relative shifts 
between the two spectra are less than 1 mA except in one case 
where it is 1.7 mA. We note that the quadratic refinement tech- 
nique (instead of a Gaussian fitting) also gives similar results. 
To derive the global trend of the relative shift, we have extended 
our cross-correlation analysis, to the entire wavelength range. 
The result of the analysis is shown in Fig. [2] The shifts are ob- 
tained in the same way as in Fig. ^ The average of the mean 
relative shifts over the entire wavelength range is 0.01 mA with 
an rms deviation of 1 .09 mA. In what follows we investigate 
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Fig. 1. The points with error-bars are the cross-correlation co- 
efficients plotted as a function of the relative shift between 
UVES and HARPS Th-Ar spectra. The cross-correlation is per- 
formed using groups of five consecutive unblended emission 
lines in the vicinity of the different metal absorption lines of 
the z a b s = 1.1508 system. The observed wavelength of the re- 
gion around the metal line is given in the top left corner in each 
panel. The dotted line is the best Gaussian fit to these coeffi- 
cients. This is used to derive the relative shift between the two 
spectra with sub-pixel accuracy. The mean relative shift and l<x 
error as well as the central wavelength of the region used are 
given in each panel. 

the absolute wavelength calibration accuracies of the two in- 
struments. 

3.2. Testing absolute wavelength calibration error of 
UVES and HARPS 

To test the absolute wavelength calibration accuracy we com- 
pare the central wavelength of strong un-blended emission lines 
in the extracted Th-Ar lamp spectrum with the wavelengths tab- 
ulated in Cuyper et al. (1998). We model the emission lines by a 
single Gaussian function. The best-fit line-centroid along with 
other parameters of the models and errors are determined by a 
X 1 minimization procedure. In many cases we find it difficult to 
fit the lines with reduced^ 2 as 1. In such cases we have scaled 
the flux errors by square root of the reduced x 1 and re-run the 
fitting procedure. In this way, we have avoided any underesti- 
mation of the errors on the best fit parameters, assuming that 
the actual errors on the flux of the Th-Ar lamp spectrum was 
somehow underestimated. 

The difference between the best-fit line centroid, in the ex- 
tracted lamp spectra and the wavelength quoted by Cuyper et 
al. (1998) is plotted in Fig.fJ] The wavelength range shown in 
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Fig. 2. The left panel shows the shift between HARPS 
and UVES lamp spectra derived by performing the cross- 
correlation (as shown in Fig. [Q over the wavelength regions 
consisting of five consecutive unblended Th-Ar lines. The his- 
togram of the mean shift is shown in the right panel. The mean 
shift is 0.01 mA and the rms <x = 1.09 mA. 
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Fig. 3. A/1, the offset of the centroid wavelength of the emission 
lines in the Th-Ar lamp spectra (obtained using Gaussian fits) 
with respect to the wavelengths given by Cuyperet al. (1998) is 
plotted versus wavelength. The left-hand side upper and lower 
panels show the results for the HARPS and UVES spectra re- 
spectively. The corresponding right-hand side panels provide 
histograms. The root mean square deviation of AA around zero 
(cr) is stated explicitly. 



this figure is the one covered by the main Fe n and Mg n lines 
of the z a bs = 1-1508 system. We find the rms of the deviation 
(AA in Fig. [3) around zero to be, respectively, 0.87 mA and 
4.08 mA for the HARPS and UVES lamp spectra. This clearly 
demonstrates that the shifts between the HARPS and the UVES 
lamp spectra measured from the cross-correlation analysis (i.e 
< 1 mA) are well within the wavelength calibration accuracy 
of UVES. 
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Fig. 4. The spectral resolution (R = A/FWHM) measured 
from Gaussian fitting of emission lines in the Th-Ar lamp spec- 
tra is plotted versus wavelength. The left-hand side upper and 
lower panels show the result from HARPS and UVES spectra 
respectively. The corresponding histograms are shown in the 
right-hand side panels. The mean along with the standard devi- 
ation are stated explicitly. 



In addition, we have used the best-fit FWHM of the 
Gaussian fit of the lamp lines to derive the spectral resolution 
(R = A/FWHM) of the spectrum. The resolution measure- 
ments are shown in Fig. |4] The mean resolution and standard 
deviation for HARPS and UVES are found to be R = 1 12, 200 
and cr = 8, 400; R = 55, 100 and cr = 7, 600 respectively. 

3.3. Effect of calibration error on Aa/a measurement 

Next we investigate how the scatter in wavelength calibration 
(AA) translates into a scatter in Aa/a . We follow the method 
used by Murphy et al. (2003) for this purpose. We randomly 
choose 3 emission lines in the lamp spectrum, with a rest wave- 
length close to each of the observed wavelengths of the Fe n and 
Mg ii lines used in the analysis of the variation of a. There are 
two Mg ii lines, ,12796 and ,12803, and five Fe n lines, ,12344, 
,12374, /12382, ,12586, and ,12600. Thus we have 21 (7x3) 
lines per realization. By choosing 3 lines, we mimic 3 distinct 
components in the actual absorption system. We assume that 
the measured shift in the emission line centroid away from the 
actual value is caused by the variation in a. To estimate this 
variation, we use the analytic fitting function given by Dzuba 
et al. (2002), 



w + qx. 



(1) 



Here, w a and w are, respectively, the vacuum wave number (in 
units of cm"') measured in the laboratory and the modified 
wave number due to a change in a; x = (Aa/a + l) 2 - 1 and 
q is the sensitivity coefficient. At each chosen lamp emission 



H. Chand et al.: Variation of the fine-structure constant 



5 



a = 


0.086x10-' (UVES Only Fe II) - 

_ t tL t if . . 1 rh i i Jf: 
















15 30 45 



(T = 0.019x 10- 5 (HARPS Mg II/Fe II) 



15 30 



0.2 




t7 = 6.096x lO" 5 (UVES Mg Il'/Fe II) ' ' ' " 


; 11 tit tYl tftf 1 


M ^AP Jt } if 4 V 








40 80 
Realisation 



15 30 
Number 



Fig. 5. Inferred Aa/a due to uncertainties in the wavelength 
calibration. The results presented in the middle and bottom 
panels are when we consider 5 main Fe n lines along with 
the Mg ii doublet in our analysis. For each realization a set 
of 21 Th-Ar lines are randomly chosen (3 in the neighbor- 
hood of each 5 Fe n and 2 Mg n lines of the z a bs =1.1508 
system). Aa/a is computed from the measured deviations in 
the line centroid by assigning the q coefficient (Dzuba et al. 
2002) of neighboring metal line to the Th-Ar emission line. 
The left-hand side middle and lower panels show the result for 
the HARPS and UVES spectra respectively. The histogram for 
both cases are shown in the right-hand side panels. The top 
panels gives the results for the UVES data when we consider 
6 Fe ii lines (i.e Fe H/11608 and 5 main Fe 11 lines) in the anal- 
ysis. The cr of the distribution refers to a typical error on the 
measurement of Aa/a due to wavelength calibration alone in a 
single system with 3 distinct components. 



line we assign the q value of the neighboring metal absorption 
transition. 

All the lamp emission lines in each realization are fitted 
simultaneously with Gaussians, for one fixed value of Aa/a . 
Here, the Aa/a value is used to modify the rest wavelength of 
the emission lines using the q coefficients given by Dzuba et 
al. (2002) for the corresponding metal lines. This procedure is 
repeated for a range of Aa/a , from -2.0 x 10~ 5 to 2.0 x 10~ 5 
in steps of 0.02 x 10~ 5 to achieve^ 2 as a function of Aa/a . 
The x 1 versus Aa/a curve is used to extract the best fitted 
Aa/a (with error-bars) in a similar way as is used in the ab- 
sorption system (discussed in the next Section). The measured 
spurious Aa/a for 100 random realizations are plotted in Fig. [5] 
both for HARPS (left-hand side middle panel) and UVES (left- 
hand side lower panel) lamp spectra. In the top panel we give 
the results for similar analysis of UVES spectrum considering 
6 Fe 11 lines (i.e including Fe 11/U6O8 instead of Mg 11 doublet) 
alone. 



We notice that the measured values of Aa/a obtained in this 
experiment have a Gaussian-shape distribution with cr of 0.02 x 
10~ 5 for HARPS and cr = 0.1 x 10~ 5 for UVES. As the sys- 
tem under consideration is known to have much more than 3 
components, the above quoted values are conservative errors 
due to uncertainties in the wavelength calibration. Murphy et 
al. (2003) have also carried out such analysis for HIRES Th-Ar 
lamp spectra. Their weighted mean from the sample of 128 sets 
of Th-Ar lines resulted in (Aa/a) ThAr = (0.4 + 0.8) x 10~ 7 . If 
one assumes a Gaussian distribution for the individual values, 
then the central limits theorem implies that the typical cr from 
one set of Th-Ar lines in the case of HIRES should be around 



0.09 xl0~ 5 (=0.8x10" 



, which is similar to our value 



for UVES Th-Ar lamp spectra (i.e cr = 0.1 x 10~ s ). 

3.4. Effect of using different Th-Ar line tables on 
wavelength calibration 

Th-Ar reference wavelengths are taken from the compilations 
of Palmer et al. (1983) for Thorium lines and Norlen et al. 
(1973) for Argon lines. The line lists built from these com- 
pilations and commonly used for echelle spectroscopy calibra- 
tion are available on the web-pages of the European Southern 
Observatory (ESO 1 ) and the National Optical Astronomy 
Observatory (NOAO 2 ). The two tables differ slightly, because 
the ESO Th-Ar line table is not accurate up to 4 decimal places 
as is the case with NOAO Th-Ar line table. For the extraction of 
UVES lamp spectra we have used the Th-Ar line table provided 
by NOAO. To investigate whether the use of ESO table could 
induce systematic shifts in Aa/a , we have also extracted the 
same UVES Th-Ar lamp spectrum using the Th-Ar line table 
provided by ESO. We fit a Gaussian function to the un-blended 
Th-Ar line as described in sub-section l3.2l and get the deviation, 
SAfit, of the best-fit centroid with respect to the corresponding 
value in the NOAO Th-Ar table. The deviation (6Af it ) is plotted 
in Fig.|6]as a function of the difference in the wavelengths tab- 
ulated by ESO and NOAO, AA, a b- If the wavelength uncertain- 
ties caused by the inaccurate wavelengths listed in ESO Th-Ar 
table for some of the Th-Ar lines are larger than the errors al- 
lowed by the dispersion solution, then we expect a correlation 
between SAfu and AA ta b- The lack of such a correlation and the 
larger scatter of SAfu compared to AA la i, in the figure, show 
that the effect of inaccurate rest-wavelengths of a few lines in 
the ESO line list is negligible. 

To complement this, we perform the cross-correlation be- 
tween the lamp spectra calibrated using the two wavelength ta- 
bles. The cross-correlation is performed in a similar way as 
described in sub-section 13. II Here we have shifted the UVES 
lamp spectrum calibrated using the ESO Th-Ar table over the 
same lamp spectrum calibrated using the NOAO Th-Ar line ta- 
ble. The result of the cross-correlation is shown in the upper 
panel of Fig. [71 From the figure it can be seen that the relative 
shift is not completely random. However the relative shift is 
most of the time less than 2mA and even 1mA, which is well 
within the UVES calibration accuracy. 



1 http://www.eso.org/instruments/uves/tools/tharatlas.html 

2 http://www.noao.edu/kpno/specatlas/thar/thar.html 
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Fig. 6. The difference, 6Afi t , between the best- fit centroid of 
Th-Ar lines seen in the UVES Th-Ar lamp spectrum and their 
wavelength listed in the NOAO Th-Ar table is plotted versus 
the difference between the corresponding wavelengths listed in 
the ESO and NOAO Th-Ar line tables (AA, ab ). The figure shows 
that (i) the scatter of AAtab is about a factor 3 smaller than that 
of 6 A fn (ii) no clear correlation is seen between AA, a b and SA/u . 
As a result the calibration errors due to differences in wave- 
lengths given in different Th-Ar tables is negligible as compare 
to the wavelength calibration accuracy of the instrument. 



We also repeat the exercise to derive how these wavelength 
calibration uncertainties translate into Aa/a as described in 
detail in sub-section 13 . 31 for the case when one uses for cali- 
bration the ESO Th-Ar line table (Fig.|5]for UVES lamp uses 
NOAO table). The result is shown in the lower left-hand side 
panel of the Fig.0for 100 realizations. The histogram shown in 
the lower right-hand side panel shows that the fiducial Aa/a is 
distributed like a Gaussian. As a result, we can conclude that 
the Aa/a measurements in the literature (Chand et al. 2004 & 
2005, Quast et al. 2004) using the ESO Th-Ar line table, should 
not be significantly affected by this possible systematic effect. 



4. Analysis 

In this section we present the results on the measurement of 
Aa/a using the HARPS and UVES spectra. The details of the 
analysis used here, validation of the procedure using simulated 
spectra and the error budget from x 2 analysis can be found in 
Chand et al. (2004, 2005). Here, we mainly concentrate on (i) 
comparing the methods used by Chand et al. (2004, 2005) to 
derive Aa/a with that used by Quast et al. (2004) and (ii) un- 
derstanding the effect of the decomposition of the absorption 
profiles into multiple narrow Voigt-profile. 



Fig. 7. The figure shows the effect of calibration using two dif- 
ferent Th-Ar line tables: (i) provided by NOAO and usually 
used in IRAF (ii) provided by ESO and used in MIDAS. The 
upper left panel shows the mean shift of Th-Ar lamp spec- 
trum calibrated using ESO Th-Ar line table with respect to 
the same Th-Ar lamp spectra but calibrated using the NOAO 
Th-Ar line table. The mean shift is derived by performing the 
cross-correlation as is shown in Fig.^over a wavelength region 
consisting of about 5 consecutive unblended Th-Ar lines. The 
histogram of the mean shift is shown in the right panel. The 
lower left panel shows the similar plot as in the lower left panel 
of Fig. |5] except that here Th-Ar lamp spectrum is calibrated 
using the ESO Th-Ar lines table rather than the NOAO Th-Ar 
lines table. The bottom right panel shows the histogram of the 
Aa/a values. 



4.1. Re-analysis of the red sub-system in the UVES 
data 

In the analysis of Chand et al (2004, 2005) Aa/a is not ex- 
plicitly used as fitting parameter. Instead^- 2 versus Aa/a curve 
is used to get the best fitted value of Aa/a . However, Quast 
et al. (2004) use the Voigt profile analysis keeping Aa/a also 
as a fitting parameter in addition to A', b and z. Chand et al. 
(2005), using analytic calculations, have shown that both the 
approaches should give the same result. Here we check this by 
re-analysing the absorption lines of the z a b s = 1.1508 system 
toward HE 0515-4414 using x 1 versus Aa/a curve. 

The absorption lines of this system is spread over about 
730 km s _1 (Quast et al. 2004). We have divided the whole sys- 
tem in two well detached blue and red sub-systems. The blue 
sub-system covers the velocity range -570 to -100 km s _1 and 
the red sub-system covers the velocity range -20 to +110 
km s _1 with respect to z a b s =1.1508. Our best fit Voigt-profiles 
to the blue and red sub-system using the UVES spectrum, is 




Fig. 8. The figure shows the velocity plot of observed profiles (data points with error-bars) together with the best fitted Voigt- 
profile for Aa/a = over plotted as a solid curve, to the blue (left-hand side panels) and red-subsystem (right-hand side panels) 
of the z a bs =1.1508 in the UVES spectrum. The dotted and dashed vertical lines are respectively the locations of the individual 
components obtained in this study and that of Quast et al. (2004). 



shown respectively in the left and right-hand side panels of 
Fig. [8] The vertical dotted lines are best fitted velocity com- 
ponents obtained in this study and the long dashed vertical 
lines mark the velocity components of the Quast et al. (2004). 
Apart form the component around ~ 90 km s _1 , we find al- 
most perfect matching between the components obtained with 
two different fitting codes. The variation of x 2 as a function 
of Aa/a using this initial fit (Fig. [8ji is shown in the left- 
hand side panel of Fig. [5] The scatter seen in the x 2 curve 
is mainly due to low column density of many components in 
blue sub-system (see the discussion in Chand et al. 2004). The 



position of the minimum in the x 2 curve remains uncertain 
till either we smooth the curve or fit some smoothing poly- 
nomial to it. Therefore we have fitted a polynomial function 
of 4th order minimizing the rms deviation. The best fit of the 
X 1 curve is shown by the solid line (left-hand side panel of 
Fig.|9jl. Its minimum gives Aa/a - (0.10 ± 0.22) x 10~ 5 , us- 
ing + 1 statistics. The derived position of the minimum 
does not depart significantly when we use a 2nd or 3rd or- 
der polynomial fit to the x 2 data points. Our best fitted value, 
Aa/a = (0.10 + 0.22) x 10" 5 , is very much consistent with that 
obtained by Quast etal. (2004) (Aa/a = [0.01 + 0.17] x 10~ 5 ). 
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Fig. 9. The dotted curve in both the panels indicates the variations of y 1 as a function of ha I a as measured using the UVES (left 
panel) and HARPS (right panel) spectra. The solid curve are the polynomial fit to these curve obtain using rms minimization to 
avoid local fluctuations. Dark rectangles with error bar indicate the position of the minimum with one sigma error-bar obtained 
fmvaxlm + 1 statistics. The^ 2 curve in the left-hand side panel is derived based on the initial fit of UVES data shown in Fig. [8] 
while the curve in the right-hand side panel is obtained by using simultaneously the initial fit of HARPS data shown in Fig. 1121 
and Fig.lT3*l 
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Fig. 10. Comparison of parameters derived for individual components from the fit of the UVES spectrum in this study and that 
of Quast et al. (2004). 



The best fitted column densities and Doppler parameters in in- 
dividual components also agree well (see Fig. llOt . The larger 
errors in the measured quantities in the present study is mainly 
due to higher values of the error assigned to the flux in individ- 
ual pixels. Thus the analysis presented here clearly shows that 



the analysis used by us in Chand et al (2004, 2005) produces 
consistent results. 

In addition we have also performed the analysis of UVES 
spectra by excluding the weaker Fe n lines from the blue sub- 
system and heavily saturated strong Fe nvi/12383,2600 lines 
from the red-subsystem, (see discussion in Chand et al. 2004). 
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Fig. 11. Absorption profiles in the red sub-system of the z abs = 1.1508 DLA toward HE 0515-4414 as observed with HARPS 
and UVES plotted on a velocity scale. The normalized UVES spectrum is shifted in the y-direction by one unit for the sake of 
clarity. The data points with error-bars correspond to the observed spectra. Over plotted as a solid curve is the best Voigt-profile 
fit based on the UVES data alone (same as in right-hand side panels of Fig. [8jl. For HARPS data the fit based on UVES data 
has been convolved with HARPS instrumental profile. The figure demonstrates the requirement for extra components, as evident 
from the higher spectral resolution HARPS spectrum (see for example region around -20 to 30 kms~ l ). 



In this case the x 2 curve is found relatively less fluctuating as 
compare to the left-hand side panel of Fig. [9] and has resulted 
in Aa/a = (0.00 ± 0.26) x 10~ 5 . 



4.2. Aa/a from the HARPS data 

The decomposition of the absorption profiles in sub- 
components is expected to be better defined from the HARPS 
spectrum because of its superior spectral resolution. In Fig. II II 
we compare the profiles of the Fe n lines in the red sub- 
system as observed with HARPS and UVES. The best multi- 
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Fig. 12. Absorption profiles in the blue sub-system of the z a b s =1.1508 on a velocity scale. The normalized UVES spectrum is 
shifted by unity for better visualization. The data points with error-bars correspond to the observed spectra. Over plotted as solid 
curves are the best Voigt-profile fits at Aa/a = 0. The Voigt-profile fits of both the HARPS and UVES data, shown here are based 
on the component structure derived by imposing the condition that the HARPS (R = 1 12000) as well as the UVES {R = 55, 000 
and better S/N) data should be fitted with parameters that are consistent with each others. The dotted vertical lines mark the 
positions of components required to fit the HARPS data. The thick ticks mark the position of components as is derived using the 
UVES data alone (Fig. |3l 



component Voigt-profiles fit using the UVES spectrum alone 
is over plotted. To fit the HARPS data we need additional 
components, as is apparent in the region around -20 to +30 
km s where consistent differences are seen for all profiles 
between the HARPS spectrum and the fit to the UVES data 
alone. However, the UVES spectrum has the advantage of hav- 
ing higher S/N. Thus, in our analysis we fitted simultaneously 
both HARPS and UVES data using the same component struc- 
ture and the appropriate instrumental functions. We initially fit- 
ted the HARPS data and used the derived parameters to fit the 
UVES data. The process was repeated until the residuals along 
the profiles are symmetrically distributed around zero and the 
best-fit parameters from these two data sets are consistent with 
one another within measurement uncertainties. In this exercise 
we have not included the line Fe n/11608 (covered only in the 
UVES spectrum) so that our derived component structure is not 
artificially bias towards Aa/a = 0. 



Our best-fit Voigt-profile components that simultaneously 
fit the HARPS and UVES spectra are shown in Fig. I12ll3l re- 
spectively for the blue and red sub-systems. The best-fit pa- 
rameters are listed in Table. ^ The component identification 
number (C.N), redshift (z), velocity dispersion (b), and Fe n 
column density (AO, for each component are listed respectively 
in columns 1, 2, 3 and 4. The last column of the table lists the 
relative velocity of the components with respect to z a b s = 1 .1508. 
We find that the blue and red sub-system (Fig. 1 1 21 1 3I > require 
respectively 3 and 6 extra components compared to the min- 
imum number required to fit the UVES spectrum alone with 

X 1 = i. 

We evaluate the best-fit Aa/a value using the high res- 
olution HARPS spectrum for the five main Fe n lines and 
the UVES spectrum for Fe H/11608 considering both the blue 
and red sub-systems simultaneously. Here it should be noted 
that the Fe ii/U608 is crucial for Aa/a measurement due to 
its opposite sensitivity for Aa/a (negative q coefficient) com- 
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Fig. 13. Same as Fig.^]for the red sub-system. In addition, the figure also illustrate more clearly that the HARPS data require 
more components (15 components, shown by dotted line) compared to the UVES data alone (9 components, shown by thick 
ticks). Also note that the fit shown for the UVES data is based on the component structure obtained in conjunction to HARPS 
data. For comparison the component required to fit the UVES data alone (from Fig.[8]right-hand side panels) are marked by thick 
tick. 



pared to the other main Fe n lines. However as its observed 
wavelength range (« 3460A) is not covered by the HARPS 
spectral coverage (3800 - 6900A), we have to use it from the 
UVES spectrum for constraining the Aa/a value. The^ 2 ver- 
sus Aa/a curve is shown in the right-hand side panel of Fig. [9] 
The scatter seen in the x 2 curve is mainly due to the low S/N 



ratio and low column density of many components as can be 
seen from Tabled (see the discussion in Chand et al. 2004). 
The continuous curve gives the 4th order polynomial fit to 
the x 1 data points using rms minimisation. Its minimum gives 
Aa/a = (0.05 ± 0.24) x 10~ 5 , using x 2 min + 1 statistics. This 
result is consistent with the Quast et al. (2004) measurement 
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(Aa/a = [0.01 ± 0.17] x 10~ 5 ) based on the UVES spectrum 
and lesser number of components. Thus in this particular case 
lack of information on the additional components in the UVES 
spectrum does not seem to affect the final result. 

5. Result and discussion 

In this paper, we present a very high resolution (R = 1 12,000) 
spectrum of QSO HE 0515-4414 obtained using HARPS. We 
have used the high wavelength calibration accuracy and high 
spectral resolution capabilities of HARPS to address the fol- 
lowing issues. 

We compare the lamp spectra obtained with UVES and 
HARPS. Using cross-correlation analysis we show that any 
possible relative shift between the two spectra are within 2 mA. 
Using Gaussian fits to unblended lamp emission lines, we find 
that the absolute wavelength calibration of HARPS is very ro- 
bust with rms deviation of 0.87 mA with respect to the wave- 
lengths tabulated in Cuyper et al. (1998). This is about a fac- 
tor of 4 better than that of UVES (cr = 4.08 mA, see Fig. 0. 
Thus the small shifts noted between the HARPS and UVES 
lamp spectra are well within the typical wavelength calibration 
accuracy of UVES. We have derived the error on Aa/a mea- 
surements due to the calibration accuracy alone. For UVES and 
HARPS spectra this is found to be respectively cr = 0.96 x 10~ 6 
and cr = 0.19 x 10~ 6 for a typical system with three well de- 
tached components. The value obtained for the UVES spectrum 
is also consistent with that of HIRES (Murphy et al. 2003). 

This shows that HARPS is the ideal instrument for this kind 
of measurement. Unfortunately it is mounted on the 3.6 m tele- 
scope at La Silla and only HE 0515-4414 is bright enough to 
be observed in a reasonable amount of time. This shows as well 
that the UVES spectra reduced (or calibrated) with the UVES 
pipeline and used in the literature to constrain Aa/a (Srianand 
et al. 2004 and Chand et al. 2004, Quast et al. 2004, Chand 
et al. 2005) do not suffer from major systematic error in the 
wavelength calibration. 

We have obtained the accurate multi-component structure 
using the higher resolution data (R * 112,000 for HARPS 
compared to » 55,000 for UVES). The best fit to the pro- 
files obtained by fitting simultaneously the HARPS data (of 
higher resolution) and the UVES data (of better S/N ratio) re- 
quire additional components as compared to the fit using the 
UVES data alone (Quast et al. 2004). Using this new sub- 
component decomposition and both HARPS and UVES data, 
we find Aa/a = (0.05 ± 0.24) x 10~ 5 . This is consistent with 
the results derived by Quast et al. (2004) from the UVES data 
alone. Indeed, we have in addition re-analyzed the UVES data 
which was used in Quast et al. (2004) (without using the com- 
ponent structure from HARPS data), to estimate the effect of 
different independent algorithms used to obtain error spectra, 
to combine the data, to fit the continuum and to fit the absorp- 
tion lines. We find that the best-fit parameters as well as the 
Aa/a measurement (Aa/a = [0.10 + 0.22] x 10~ 5 ), obtained 
by our independent analysis, are consistent with that of Quast 
et al. (2004) (Aa/a = [0.01 ± 0.17] x 10~ 5 ). 

We note that the precision on the Aa/a measurement ob- 
tained using the HARPS spectrum, which is of high resolution 



and low S/N ratio, is similar to that obtained from the UVES 
spectrum, which is of lower resolution and higher S/N ratio. 
Therefore, the improvement in the wavelength calibration ac- 
curacy by an order of magnitude using HARPS will be effec- 
tive to improve the constrain on Aa/a only if high S/N ratio 
can also be obtained. This could be possible if an instrument 
such as HARPS can be mounted on bigger telescopes. 
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Table 1. Results of the Voigt profile fit of Fe n lines at z a bs= 1.1508 toward HE 05 15-4414. 
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' relative velocity with respect to z a b s = 1.1508. 

' The redshift (z) of these components are kept fixed. 

' The Doppler parameter, b, of these components are kept fixed. 



